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Abstract 

Suppose the lightest superpartner (LSP) is observed at colliders, and WIMPs 
are detected in explicit experiments. We point out that one cannot immedi- 
ately conclude that cold dark matter (CDM) of the universe has been ob- 
served, and we determine what measurements are necessary before such a 
conclusion is meaningful. We discuss the analogous situation for neutrinos 
and axions; in the axion case we have not found a way to conclude axions are 
the CDM even if axions are detected. 
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I. INTRODUCTION 



Despite the many successes of the standard cosmology, we still do not know the compo- 
sition and the amount of energy density in our universe (for a review, see for example J^j). 
The existence of matter which does not emit much light is certain from the fact that stars 
and other luminous matter contribute only a tiny fraction (about 0.004) || of the critical 
energy density (required for a spatially flat universe) while the orbits of stars around galaxies 
(see e.g. [Q) indicate that the gravitating energy density is about ten times larger .Q From a 
cosmological point of view, a variety of observations, such as the peculiar motions of galaxies 
as well as the masses of clusters of galaxies, corroborate the existence of dark energy and 
indicate that Qm, the ratio of the average matter-energy density in our Hubble volume to 
the critical energy density, is about 1/3 if the Hubble expansion rate is about 70 km/s/Mpc 
as the current observations indicate (see e.g. j|]). Although without other constraints this 
dark energy density could be in the form of ordinary baryonic matter (objects made of neu- 
trons and protons), the standard model of big bang nucleosynthesis (BBN) calculations and 
the measurement of primeval abundance of deuterium gives strong evidence that the baryon 
density is VLbIi 2 = (0.02 ±0.002) which implies Qb of about 0.04. Hence, most of the matter 
energy density in the universe seems to be in the form of nonbaryonic dark matter (NBDM). 

Measurements of cosmological parameters will determine ^Inbdm accurately, but are 
unlikely to help us know the actual character of the NBDM, since they mainly depend on its 
gravitational interactions. Some information on NBDM collisional energy loss and clustering 
will also come from astrophysics information. For perhaps the most likely forms of NBDM 
we can hope to observe the actual particles in laboratory experiments, and calculate to a 
few percent accuracy the actual contribution to Q ||. 

Remarkably, the more likely extensions to the standard model of particle physics provide 
candidates for the nonbaryonic dark matter, candidates which existed even before the need 
for NBDM was established. Indeed, given that these extensions to the standard model are 
theoretically compelling from particle physics considerations, without cosmological consid- 
erations, their provision of dark matter energy of the right order of magnitude abundance 
provides an independent hint at the existence of physics beyond the standard model. Among 
the various candidates for this nonbaryonic dark matter (see e.g. ) , a typical unified super- 
symmetric theory with conserved R-parity will have a stable lightest supersymmetric particle 
(LSP) that will probably constitute most of the cold dark matter (CDM). In addition, the 
non-baryonic dark energy would most likely consist of neutrino hot dark matter (z/HDM), 
axionic CDM (ACDM), and the cosmological constant (A). Given a Lagrangian for the 
unified theory (UT) and thermal equilibrium initial conditions for the fluid determining the 



In this paper, we will use the standard notation of ratio of energy density in X to the critical 
energy density as fix ■ Since the critical energy density is determined by the Hubble expansion rate 
which is uncertain, we will use the usual parameterization H = MOO km/s/Mpc and often write 
the energy density as Q.xh 2 value instead. 
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cosmological evolution,^ one can calculate the dark matter content of our universe by solving 
the Boltzmann equations. 

Many other NBDM candidates have been proposed ||. In this paper we will only analyze 
the situations for neutrinos, axions, and LSPs in some detail. Similar conclusions hold for 
the rest. At appropriate places in the paper we will consider non-thermal evolution. 

The neutrino situation is simpler. Once the cosmology is known and the neutrino masses 
are determined, one can compute the relic number of neutrinos, multiply by the masses, and 
obtain Vt u . If neutrinos have anomalous interactions they have to be included. 

The axion case is difficult. Assume axions are observed in terrestrial detectors, and 
perhaps even information about them comes from astrophysical data. Then the thermal 
density Q a (t) can be computed precisely enough @. However, coherent oscillations of axion 
zero modes give non-thermal contributions fl a (nt) that depend strongly on a "misalignment" 
factor, related to the value of the axion field at the confinement transition. We are aware of 
no way to determine fl a (nt) , which can be the dominant contribution, so we conclude that it 
may not be possible to even know how much of the CDM is £l axion = + Q a ( n t)- 

The LSP case turns out to be difficult, but solvable. We first demonstrate that knowing 
the LSP mass, and even knowing its cross section on a nucleon, do not allow one to know 
its contribution to Qlsp to better accuracy than an order of magnitude or so. Part of this 
uncertainty arises from a current lack of knowledge of the phases of the soft supersymmetry 
breaking Lagrangian (C so f t ), but even if one arbitrarily took all phases to be zero or n at 
least a factor of six uncertainty in Qlsp remains. And for an important question such as the 
composition of the CDM, one would not want to make unwarranted assumptions about the 
phases or other relevant parameters. We then show that knowing the relevant parameters 
of Csoft, and tan/3, to about 5% allows a determination of Qlsp to similar accuracy. Such 
measurements will be possible by combining data from hadron colliders, low energy (e.g. 
electric dipole moments) experiments, B-factories, and a lepton collider with a polarized 
beam and sufficient energy to produce several superpartners (and appropriate luminosity). 

If WIMPs are observed in explicit (i.e., non-collider) detection experiments (direct un- 
derground WIMP scattering on nuclei, or space-based, or indirectly via v interactions from 
WIMP annihilation in the sun or earth) there are additional large uncertainties which we 
do not consider in this paper. These include going from WIMP-nucleus to WIMP-quark 
cross sections, local density and velocity distributions of WIMPs, how long antiprotons or 
positrons can persist in the galaxy, and so on. All of our results hold even if these other 
factors could be controlled. 

In the next two sections we discuss the LSP case in detail, first analytically and then 
numerically. Then we turn to some cosmological considerations, and examine the axion 
case. In all cases we assume that h 2 will be known accurately before it is needed to achieve 
a particle physics knowledge of Qnbdm, so we do not include errors in h 2 . We also assume 



2 The cosmological initial conditions can be determined by an inflationary model, which in principle 
can also be determined by the UT if the initial conditions for inflation are set by an even more 
fundamental principle or UT offers a unique inflationary history. We assume that the underlying 
cosmology is known when the relic abundances are finally calculated. 
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that the usual Boltzmann treatment is sufficient || to estimate the possible uncertainties 
and that refinements of the calculational procedure necessary for more accurate calculation 
can be accomplished when appropriate. For examples of such refinements, see Ref. || and 
references therein. 



II. ANALYTIC ESTIMATES 



In this section, we give an analytic estimate for the uncertainties that can be expected 
in the calculation of LSP CDM, given that some of the parameters in the MSSM has been 
measured to a certain accuracy. The most common situation for most parts of the MSSM 



parameter space is when no other particle mass is within m x /20 of the LSP mass m x [10,1] 



In that case, self-annihilation determines the relic abundance. Assuming as usual that 
the annihilation products are in chemical thermal equilibrium with the massless degrees of 
freedom in a radiation dominated flat FRW universe, the simplified Boltzmann equation 
governing the relic abundance can be written as 



dx 



fl 



(1) 



where / = n x /T 3 is the LSP volume density scaled by the cube of the temperature which 
sets the scale for the volume density of the photons, (av) is the thermal averaged self 
annihilation cross section, x = T/m x is the scaled temperature, g* counts the degrees of 



freedom contributing to the entropy, and f 



-3/2 



27T 



-1/as ; 



is the nonrelativistic approximation 



of the thermal equilibrium volume density of LSPs. Starting from a thermal equilibrium 
initial conditions, / tracks f until the freeze-out temeperature Xp and then / decouples 
from f . Approximating this decoupling to occur sharply at a particular temperature xf , 
one finds an expression for the relic density as 



n 



r (/ {av)dx) lp c 



(2) 



pi Jx o 



where p c is the critical energy density that only depends on the cosmology. The approximate 
expression for the freezeout temperature can be written as 



1 



Xp 



]n[m x £{av)] + \ \tlxf 



(2tt)3\ 2g.. 



Propagating the error in quadratures, the fractional error is characterized by 
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where 5P{ denotes the uncertainty in parameter Pj. Neglecting the T and (^uncertainties, 
we can write 



A,- = 



Ixo y/<h((Tv)dx 



where we have used the fact that sj? < 1. We can approximate d ^ ~ (av) -p; where 
Ti ~ 0(1) (which is a good approximation in the region of the parameter space where the 
Pi dependence is analytic) and conclude 

A,; ~ — . 



Hence, we conclude that 



5P 



where the strength of the errror contribution is determined by = ^j^p^ which is what we 
expect. 

For example, consider a typical nonresonant self-annihilation cross section of an LSP 
going to two fermions through a sfermion exchange. We have 



(av) 



„ g2 ffx 



64tt\ 
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m. 



mi 



Ciirij + C2xm 2 x 



{m 2 - + w? x — m 2 ) 2 



f 



where C\ and C2 are 0(1) dimensionless constants and g 2 j x is the fermion-sfermion-neutralino 
coupling. Then considering the contribution of m x to the fractional error, we find to leading 

■i . rat rat 

order m — ~ j — 



ra x raj 



2 1- 



l + (^) 2 



which is about < 1 . Considering the contribution of mj to the uncertainty, we find similarly 
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1+ iHbc 

which has a magnitude of about > 2. Finally we see that the coupling will also contribute 

T 9ffx = 2 ' 

showing that the abundance is most sensitive to the sfermion mass and the LSP-sfermion- 
fermion coupling. In addition we see that relative uncertainties in parameters are likely to 
lead to even larger relative uncertainties in Qlsp- 
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The standard scalar neutralino cross section on protons [O is expressed as 



scalar 



Am 2 T ,2 

J V 



71 



(3) 



where the effective scalar interaction coupling f p is usually dominated by the CP-even Higgs 
parton level exchange between quarks and neutralinos |13 



scalar 



jp 



m r . 



f H 2 f* 



q=u,d,s 



q=c,b,t 1 



(4) 



The matrix element coefficients and f^ G = 1 — J2 q fr q can be extracted from pion-nucleon 
scattering using chiral perturbation theory and are subject to large uncertainties which are 
reflected in the neutralino-proton scattering cross section |TJj . 



The parton level couplings, on the other hand, depend entirely on the SUSY Lagrangian 
parameters and Higgs masses 



J a 



m. 



E 



Ji)Ji) 
2 

mi 



i=l,2 '"'Hi 

with the (in general complex) neutralino couplings to the light CP-even Higgs 



.(i) 



- (gN* 12 - g'N* n ) (N* 3 sin a + N* u cos a) 



and to the heavy Higgs 



.(2) 



2 (^^1*2 - g'Nti) sin a - N* 3 cos a) , 



where a is the Higgs mixing angle. The quark-Higgs couplings depend on weak isospin 
quantum number and we have for the up type quarks 



.(i) 



g cos a 
2mw sin (3 ' 



and for the down type quarks 



g sin a 

2mw cos (3 ' 



c (2) 



g sin a 

2mw sin (3 



g cos a 
2mw cos (3 



(5) 



(6) 



(2) 

In the large limit a ~ (3 — n/2 and the c d coupling is enhanced by a factor of tan/3. 
As a result, for small values of tan/3 (less than 4) the cross section is dominated by the 
light Higgs exchange while for large values of tan (3 the heavy Higgs exchange contribution 
prevails. 

In analogy to the analysis of the relic density we can determine accuracy of the cross 
section calculation as a function of the variation in individual parameters 



(I 
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where Sj is the variation coeficient corresponding to the uncertainty 8 Pi in parameter Pi. 
Assuming that there is indeed a dominant contribution from one of the Higgs bosons, we 
find that the Higgs and quark coupling both contribute to the uncertainty in a similar way 

s Cq ^s Cx ~2 (7) 

and the Higgs mass contribution to the uncertainty is 

s mH ~ -4. (8) 

while the cross section is largely insensitive to variations of the neutralino mass since the 
reduced mass of the system is very close to the mass of the proton. 

It is important to realize that the neutralino-Higgs couplings depend crucially on the 
neutralino mixing matrix. Both gaugino and Higgsino components are required to partic- 
ipate if the couplings are not to vanish. The neutralino mass matrix depends on complex 
quantities M 1; M 2 and fi with potentially large phases |]ISJ which can significantly modify 
the lightest neutralino composition and subsequently its couplings to the Higgs bosons. It 
is important to include these complex phases in the general analysis of the relic density and 
the neutralino proton cross section in order to be able to estimate the overall uncertainty 
in these quantities which can be achieved once the SUSY Lagrangian parameters including 
the phases are measured. 



NUMERICAL RESULTS 

We illustrate the general behavior of the neutralino relic density and proton elastic 
cross section on a characteristic set of the MSSM parameters which can provide neutralino 
relic abundance in the cosmologically relevant region and a cross section small enough to 
be allowed by direct detection experiments. Since the neutralino scattering cross section 
grows with tan/3 models with low and moderate values of tan/3 can easily satisfy direct 
detection constraints and still provide a significant neutralino abundance. In all of our 
numerical calculations we use the following reference set of parameters - Mi = 80 GeV, 
m A = 250 GeV, tan/3 = 3, <p 1 = <p IM = 0, m p = 110 GeV, m ih = 125 GeV, m~ lR = 110 GeV, 
rriq L = 420 GeV, m^ R = 400 GeV, mg R = 380 GeV. 

This set has been chosen so that the resulting neutralino is predominantly a bino and the 
effects of the light Higgs and Z pole neutralino annihilation pole as well as any co-annihilation 
effects are minimized. Our choice of parameters leads to the values of Qh 2 ~ 0.148 and 
cr p ~ 11.4 x 10~ 9 pb. Since we are working in a general parametric framework, the soft Higgs 
masses can be chosen so that electroweak symmetry breaking conditions are satisfied. These 
numbers are not special, and only illustrate typical results. 

First let us turn to the discussion of the neutralino observables sensitivity to the CP- 
conserving parameters appearing in the supersymmetric Lagrangian. Obviously, the most 
important ones are Mi, M 2 , fi and tan/3 which enter into neutralino mass matrix and 
determine both the mass of the lightest neutralino m x and its Higgs and lepton-slepton 
couplings. Another significant variation comes from the change in the mass of the CP-odd 
Higgs boson which influences the neutralino-Higgs couplings and the heavy CP-even Higgs 
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FIG. 1. Variation of the neutralino relic density (a) and of the neutralino-proton elastic scat- 
tering cross section (b) with 5 = Act/ci for the most relevant SUSY parameters. 

boson mass. The dominant slepton exchange contribution to the neutralino annihilation 
cross section depends crucially on the slepton mass while the scattering cross section is 
insensitive to the scalar masses. 

Fig. 1 shows the dependence of the neutralino relic density and proton cross section 
on the relative variation of the important parameters within a ±20% range. As shown in 
frame a), the relic density is mostly sensitive to variations in the slepton masses appearing 
in the dominant annihilation diagram - as the slepton mass increases the annihilation cross 
section is more suppressed and Qh 2 increases. All the other parameters enter the neutralino 
sector and their variations are reflected in the variations of the LSP mass and the LSP- 
slepton-lepton coupling. The most significant is the variation of the density with the bino 
mass Mi which determines m x and subsequently the dominant s-wave contribution to the 
neutralino annihilation cross section. It is relatively stable with respect to changes in M2, /i 
and tan j3 since the annihilation cross section depends weakly on the neutralino composition. 
In summary, the overall variation of flh 2 is at most ±25% for any single parameter in the 
given range of the input SUSY parameters. 
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On the other hand, the spin independent part of the neutralino-proton scattering cross 
section depends crucially on the gaugino-Higgsino composition of the lightest neutralino 
since both components take part in the neutralino-Higgs interaction. This fact is reflected 
in frame b) which shows a comparable sensitivity to M%, M 2 and fi. It is obvious that for 
the given range of 5 = Ac/c the change in a p can be as big as a factor of two. Note that 
in our particular case (tan/3 = 3) the cross section is dominated by the light Higgs boson 
exchange and consequently the sensitivity to rriA is limited. As tan j3 increases, the heavy 
Higgs exchange takes over the cross section and the sensitivity to m^, which determines the 
heavy Higgs boson mass, increases correspondingly. 
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FIG. 2. Plot of regions in the £lh 2 -a p plane resulting from varying all relevant parameters from 
Fig. 1 around their standard value (square) within 5% (green (grey) region) and 20% (dots). The 
CP-violating phases are set to zero. 

In order to estimate the accuracy with which we can determine the relic density and the 
proton scattering cross section once we have measured the SUSY parameters at a collider 
experiment we plot the range of both quantities when the SUSY parameters are all varied 
within a 20 % range from the central value. Fig. 2 shows the resulting region in the Qh 2 -a p 
plane. In part of parameter space the variations of Fig. 1 can combine, and lead to much 
larger ranges of variation in Qh 2 . In particular, the points in the upper left region with 
small relic density result from non-linear effects associated with the presence of an s-channel 
light Higgs exchange annihilation pole, which can drastically increase the annihilation cross 
section, and parameter points on the tail of the Breit-Wigner resonance have small neutralino 
relic density. This is a fairly generic situation since the allowed range for neutralino relic 
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abundance imposes a limit on the neutralino mass. Note that the 5% error region allows a 
reasonable determination of Qh 2 . To understand the weak constraints from direct detection 
experiments, imagine a horizontal line across Fig. 2, with a height uncertainty coming from 
the nuclear physics and astrophysics ambiguities in extracting a p , and notice the resulting 
uncertainties in Qh 2 . 




FIG. 3. Variation of the relic abundance Qh 2 (a) and LSP proton cross section a p (b) (in 
1CP 6 pb) with CP-violating phases affecting the neutralino mass matrix. All other parameters are 
set to the standard set values. 

Next we examine the effects of soft SUSY phases on these questions. It has been demon- 
strated recently |TJj that all soft SUSY phases can be large. Even if it turned out that some 



are smaller, e.g. from future electric dipole moment experiments, others relevant here may 



not be. Also, phases have two effects - not only do they directly enter the calculations ITS 
of Qh 2 and a p , they also make measurement of other parameters such as tan /3 more difficult, 
and impossible at LEP or hadron colliders ]T7| . 

We work in a phase parametrization consistent with [TjJ and choose <f2 — ArgiMz) = 0. 
In Fig. 3 we show the range of flh 2 (frame a)) and a p frame (b)) as the two relevant phases 
entering the neutralino mass matrix (pi and ip^ are varied in their full range while all other 
parameters are kept constant. It turns out that the relic abundance has two local minima 
and maxima and the range is bigger than an order of magnitude between the minimum 
and maximum. The neutralino proton cross section also varies by more than an order of 
magnitude but it is monotonous in the ipi plane. From Fig. 3 it is clear that by 
neglecting the phases one can be missing a crucial part of the particle physics information 
needed to determine the LSP abundance. 
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COSMOLOGICAL UNCERTAINTIES AND AXIONS 



Although the focus of our paper is examining the uncertainties in our knowledge of 
dark matter coming from particle physics, here we will briefly comment on some of the 
cosmological uncertainties related to the dark matter determination. Specifically, we would 
like to emphasize a point that can easily be overlooked: only with a firm determination of 
the cosmological history can the particle physics data weigh the universe. To make this point 
concrete, let us consider some possible difficulties with the usual approach to the dark matter 
calculation which we have assumed in our paper. First we will consider the effect of energy 
density behaving like a cosmological constant today. Then, we will consider the situation 
when the dominant contribution to zero pressure dark energy is not a thermal relic. In such 
"nonthermal" cases, particle physics data generically cannot determine the "matter" energy 
density contribution because "history" or boundary conditions of the particle physics may 
not be determined by a fundamental principle or dynamics. We will choose axionic dark 
matter to illustrate this point. Such uncertainties are relevant to all forms of dark matter. 
There can also be extra entropy production which has to be included. 

Let us first consider how the existence of energy density with a negative pressure 
|18| (which we will for simplicity call p\ ) affects the relic LSP density calculations. First, 
consider the Boltzmann equations. Suppose the spatial curvature is negligible as will be the 
case after inflation. Then, we can start with the same Boltzmann equation 

+ 3Hn = -(av)(n 2 - n 2 ) 
Since conservation of entropy is still valid, we will have 

!-HT 3 (/ 2 -/o 2 ) 



where / = -M^. In the usual scenario, we assume radiation dominance and use 

1 

t= 2H 

This is still a good approximation because the temperature at which the relic density froze 
out is O(GeV) ^nucleosynthesis an ^ we rec L u i re that p^ not destroy one of the pillars of 
standard cosmology and that Q\ be a monotonic function of time. Hence, the rest of the 
formalism follows as usual for the determination of the freezeout temperature. As far as the 
evolution of the relic density after freezeout is concerned, there is no change from before again 
because the nucleosynthesis temperature is much lower than the freezeout temperature. 

Of course, if we do not require monotonic variation of p^ then the above relation between 
the time and expansion rate may no longer valid and the relic abundance calculation may 
have to be modified. However, this is unlikely because to preserve the successes of BBN, 
the models of pa then must be fine tuned such that its energy is significant at the time 



of dark matter freezeout while being negligible at the time of BBN [JjiJ. In the unlikely 
event that such tuned history is that of our universe, then the dark matter energy density 
implied by the particle physics measurements may be significantly different from what we 
have calculated in this paper. 
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Axion scenarios offer an elegant solution to the strong CP problem. Furthermore, pseu- 
doscalars resembling axions are generically expected [^] from effective theories arising from 
compactifications of string theories.^ Unlike LSPs, the axions are expected to naturally 
contribute significantly to the cosmological energy density even when the thermal relic 
abundance does not contribute significantly Even with the assumption of inflation and 
sufficiently low reheating temperature |2l] rendering possible axionic topological defect re- 



lated contributions irrelevant, one must account for the fact that generically there may be 
coherent oscillations of the axion zero modes (commonly called misalignment contribution 
0) giving a nonthermal contribution to the cosmological energy density with the magnitude 

n a(nt) h 2 = 0.13 x io ±0 - 4 ( Aqcd y^fmAeK-^-)- 1 - 18 

a(nt > V 200MeV ; J V ; lV 10" 5 eV ; 

where m a is the axion mass, / is some known monotonically increasing function accounting 
for anharmonic effects (/(0) = 1), and #i is the "amplitude" of the oscillations of the axion 
field. The value of 9\ is essentially equal to the value of the axion field at the confinement 
transition. Because there is no direct way to measure 9\ , even with a precise determination 
of m a , there will be a great uncertainty in the axionic contribution. This is in contrast 
with the most likely smaller |7j thermal relic contribution of the axion which can be quite 
precisely determined once the mass of the axion is measured. Indeed, most of the relic axion 
detection experiments such as electromagnetic spectrum observation of nearby clusters, as 
well as Sikivie-type resonant microwave cavities immersed in strong magnetic fields, mostly 
measure the axion mass, and there is no obvious direct handle on the "misalignment" angle 
Q\. Indeed, 9\ may be most likely determined randomly from the space of all possible values, 
since above the confinement transition the axions have essentially zero mass. Hence, in such 
"nonthermal" cosmological boundary condition dependent situations, the determination of 
the dark matter from particle physics data seems impossible. 



CONCLUSION 

Understanding the composition of our universe is one of the most fundamental issues we 
can study. We demonstrate that we can only be fully confident we have learned the answer 
if we can calculate from first principles and laboratory data for each type of matter X what 
Qx is, and find agreement with cosmological measurement of £l ma tter- F° r neutrinos this can 
be done if their masses can be measured, or shown to be small compared to ~ 1 eV. For 
axions we have emphasized there does not seem to be a way to do this. 

For the LSP it is possible to calculate Qlsp if measurements of soft supersymmetry 
phases, and tan/? are made. Unless such measurements are done, it is simply incorrect to 
suggest that detection of the LSP at a collider, or in a direct or indirect WIMP experiment, 
means that the cold dark matter has been observed. Indeed, there is a loose correlation of 
easier detection with smaller relic density. All measurements, of course, help constrain the 



3 However, most embeddings of axions in string compactification models run into trouble with 



scales although there may be some ways to circumvent the problem [2C|. 
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situation, and the detection of the LSP will mean that some of the CDM has been observed 
- but perhaps only a fraction. 
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